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We summarize recent work in which the shell effect, which causes the onset of 
the mass asymmetry in nuclear fission, could be explained semiclassically in the 
framework of the periodic orbit theory. We also present new results for the inclusion 
of a spin-orbit interaction in the semiclassical calculation of the level density. 



1 Introduction 

The explanation of the asymmetry in the fragment mass distributions of many 
actinide nuclei has long been a stumbleblock in fission theory. Shortly after 
the advent of the shell model, Lise Meitneru related the mass asymmetry to 
shell structure, in particular to magic nucleon numbers in one of the fragments. 
Indeed, since the classical liquid drop model (LDM) favours symmetric shapes,0 
the mass asymmetry had to be a quantum shell effect. JohanssonB attempted 
to calculate fission barriers by summing the single-particle energies of a Nilsson- 
type shell model including left-right asymmetric deformations. However, due to 
their wrong average behaviour at large deformations, he could not obtain finite 
barriers. Only_after the introduction of the ingenious shell- correction method 
by Strutinskya this became possible_Jaunching large-scale investigations of 
the shell structure in fission barriers .fflu The resulting picture of the douhle- 
humped fission barrier was nicely confirmed in many detailed experiments. u 

One interesting outcome was that the onset of the left-right asymmetry of 
the nuclear shapes starts quite early during the fission process: in nuclei like 
240 Pu, already the outer fission barrier is unstable against octupole-type defor- 
mationsrl This quantum shell effect could be related to a few specific diabatic 
single-particle states which were particularly sensitive to the asymmetric de- 
formations.E3 Today, we know thaL symmetric and asymmetric fission modes 
may coexist in the same nucleus.El Of course, the mass distributions of the 
fission fragments can only be predicted from a dynamical theory involving in- 
ertial parameters.E3 But interestingly enough, the most probable experimental 
mass ratios were found to be roughly equal to those, of the nascent fragments 
obtained statically at the asymmetric outer barrier .u t3 The onset of the mass 
asymmetry in nuclear fission was thus established as a quantum shell effect 
due to specific single-particle states in the deformed average potentials of the 
nucleons, which could not be understood in the classical liquid drop model. 
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However, new developments have taken place over the last 25 years to de- 
scribe quantum shell effects semiclassically in the framework of the periodic 
orbit theory (POT). Through the so-called trace formula, the oscillating part 
Sg (E) of the level density of a quantum system can be related to properties 
of the periodic orbits of the corresponding classical system. This approach 
was systematically developed by Gutzwiller;E£l an alternative derivation valid 
for billiard systems (i.e., particles enclosed in a cavity with_kLeally reflect- 
ing walls) was given by Balian and BlochJl3 Strutinsky et aZ.E£lO generalized 
the Gutzwiller theory and used it for the semiclassical calculation of shell- 
correction energies in various shell-model potentials (neglecting, however, the 
spin-orbit interaction). In particular, they successfully explained the system- 
atics of nucleap-ground-state deformations in terms of the leading shortest 
periodic orbits. Ii3 Similarly, the onset of the fission mass asymmetry could re- 
cently be explained semiclassicallyJHa For a comprehensive presentation of the 
POT, the trace formula and its further refinements, and its applications to 
finite fermion systems, we refer to a recent monograph.113 

In Sect. ^ we summarize the recent work on the fission asymmetry done 
by one,-Qf— lis_(M.B.) in collaboration with P. Meier, S. M. Reimann, and M. 
Sieber.EScElO In Sect. [| we present some new semiclassical results for the level 
density of a three-dimensional deformed potential, including the spin-orbit 
interaction which so far has been omitted in applications of the POT to the 
calculation of shell structure in nuclear physics. 



2 Semiclassical explanation of the onset of mass asymmetry 

In our recent semiclassical investigation of the fission mass asymmetry, Si we 
have approximated the nuclear mean field by a cavity with constant volume and 
reflecting walls, neglecting spin-orbit, Coulomb, and pairing interactions and 
considering only one kind of particles. We use the parameterization (c, h, a) 
of Ref.Q to define the boundary of the cavity. Here 2c is the length of the 
nucleus along the symmetry axis in units of the radius Rq of the spherical 
cavity (given by c=l, h=a=0), h is the neck parameter, and a ^ yields 
left-right asymmetrio_shapes. In 240 Pu, e.g., the maximum of the symmetric 
outer barrier is founcE near c=1.53, h — a = 0. nr-i 
The semiclassical trace formula for the shell-correction energy SE readal3ll3 

1 7rl 
-S po (E F ) - a po - . (1) 

The sum in (|lj) is taken over all periodic orbits (po) , but the gross-shell effects 
are dominated by the shortest periodic orbits of the system. S po = § o p • dq 



5E^Y. A po{ e f){^-) cos 



2 



© 



® 



® 




1.3 



-10 



o5E 



10 



Figure 1: Right: [Spmiclassical deformation energy SE(c, a) (h=0) in a perspective 3D 
plot (from Ref.E3). The arrows 'symmetric' and 'asymmetric' show two alternative 
fission paths. Left: shapes along the asymmetric path. Planes of the leading periodic 
orbits are shown by vertical lines (solid for stable and dashed for unstable orbits) . 

are the action integrals, T po = dS po /dE the periods, and a po are phases related 
to the number of conjugate (or focal) points along the periodic orbits. The 
amplitude A po of an orbit depends on its stability and its degeneracy; it in- 
cludes a factor that exponentially damps contributions from orbits with large 
periods T po . All quantities in (Q) are evaluated at the Fermi energy Ep. 

In Fig. [I] we show a perspective 3D plot of the deformation energy SE(c, a) 
(h = 0) obtained from Eq. (Q). The missing effect of the spin-orbit interaction 
was compensated by adjusting the only parameter of the model, Ep, such 
that the isomer minimum is found at c — 1.42, h = a = 0, corresponding to 
the quantum result for 240 Pu in the realistic Strutinsky calculations!) We see 
that the semiclassical model correctly yields the instability of the outer barrier 
against the asymmetry, and a realistic value (~ 0.13) of the parameter a at the 
asymmetric saddle point. Inclusion of the LDM energy, which varies very little 
in the deformation region shown here, does not change this result qualitatively. 
To the left in Fig. |l| we show the nuclear shapes at three points (A,B,C) along 
the adiabatic fission path. The vertical lines indicate the planes perpendicular 
to the symmetry axis in which the shortest periodic orbits are found. These are 
simply the polygons with p reflections inscribed in the circular cross section of 
the cavity with these planes. A fast convergence of SE with increasing p, i.e., 
with the length of the periodic orbits, was found in Ref.E] (p — 2 and 3 were 
essentially sufficient); other periodic orbits were also shown to be negligible. 

We thus obtain a simple physical explanation of the onset of the mass 
asymmetry of fission in terms of very few classical periodic orbits. The valley 
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of steepest descent through the deformation energy surface is simply give n0 
by the requirement that the action o£_the leading periodic orbits be station- 
ary: SS po — 0. We have also showrO that the diabatic states which cause 
this asymmetry effect quantum-mechanically have their probability maxima 
exactly on the planes containing the leading periodic orbits, and that a semi- 
classical quantization of the quasi-regular motion near these planes reproduces 
the energy levels of the diabatic quantum states surprisingly well. 

Our semiclassical picture is only qualitative. It is not meant to replace 
the Strutinsky (or Hartree-Fock) type calculations, but to give an intuitive 
physical understanding of the mass asymmetry using classical mechanics. Of 
course, one has to question if the neglect of the spin-orbit interaction can be 
justified since it is known to crucially affect the shell structure. Locally, its 
effect could here be simulated (like in RefJl3) by adjusting the Fermi energy. 
But a semiclassical calculation including spin-orbit effects is certainly desirable. 



3 Inclusion of spin-orbit interaction in POT 

We present now some new results including a spin-orbit interaction in the 
semiclassical trace formula. The cavity model of the previous section leads 
to problems with the standard (Thomas) form of the spin-orbit interaction 
which contains the gradient of the local potential. We therefore use for V(r) 
a 3-dimensional deformed harmonic-oscillator potential and add a generalized 
spin-orbit term which in the spherical limit is proportional to 1 • <x: 

H=^-p 2 + V(r) + K B(r,p)- ( T (2) 

Ira 

with 

V(r) = j Yl W M> B(r,p) = W(r)xp. (3) 

i=x,y,z 

In Eq. (B), <x = (a X: <r y , a z ) are the Pauli matrices. We allow for different 
frequencies Ui to mimic the shell-model potential of a deformed light nucleus. 
The quantum-mechanical eigenvalues of (^) were obtained by an exact diago- 
nalization in the basis of V(r) up to a maximum energy of ~ 150 ouj x . 

A first clue to the classical behaviour of this system is gained from a Fourier 
transform of the oscillating part Sg (E) of the level density from the energy to 
the time domain. To emphasize the gross-shell structure, the level density is 
Gaussian averaged over an energy range 7 < Tiuji. Due to the energy scaling 
behaviourtJ of the Hamiltonian (||), the peaks in the Fourier spectrum lie at 
the time periods T po of the classical orbits that dominate Sg (E). In Fig. || we 
show the Fourier spectra obtained for a case with irrational frequency ratios, 
both without and with spin-orbit interaction. For k = 0, we clearly see three 
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Figure 2: Squared Fourier amplitudes (in arbitrary units) of Sg (E) for the anisotropic 
harmonic oscillator (^) with frequencies lj x = 1, u> y = 1.12128, lj z = 1.25727, both 
without (left) and with spin-orbit interaction (right). Energy eigenvalues up to 
~ 50 TiLo x were included; the Gaussian averaging range was 7 = 0.2 fiu) x . 

main peaks which correspond exactly to the periods Ti = 2w/uJi of the only 
three-oeriodic orbits of this system, which are the librations along the principal 
axesO At twice these periods, three small peaks corresponding to their second 
harmonics are visible. For k = 0.1, a more complicated peak structure is 
observed, which will be explained in the following. 

The difficulty in the semiclassical description is that the Pauli operators 
in act on a two-component spinor and couple its components in a non- 
trivial way. Littlejohn and FlynnEia have developed a theory whose basic idea 
is to diagonalize the Hamiltonian (2x2) matrix locally in phase space after its 
Wigncr transformation, and to expand both Hw(?,p) and the diagonalizing 
unitary matrix in powers of h. Two problems arise hereby. First, on those 
surfaces in phase space where B(r,p) = 0, the eigenvalues become degenerate 
and the method breaks down (so-called mode conversion) . Second, one obtains 
gauge-transformation dependencies of the diagonal matrix. This problem_can 
in principle be solved by the introduction of non-caniinical coordinateso A 
more heuristic solution, proposed by Frisk and GuhrjHj leads to the determi- 
nation of the dynamics to lowest order in h by the two classical Hamiltonians 

H±(r lP ) = ^-p 2 + V(r)±K\B(r,p)\ (4) 

whose combined periodic orbits are to be_used in the trace formula. The first- 
order h corrections to the energy, in Ref .E3 denoted by \% erry and , are to 
be included perturbatively, i.e., only their contributions to the actions have to 
be included in the trace formula. This ad hoc prescription has recently been 
justified through a relatiyistic trace formula derived from the Dirac equation 
by Bolte and Keppeler.Ej From the non-relativistic reduction of their result, 
the use of the Hamiltonians (|j) and the above treatment of X^ erry and 
can, indeed, be obtained in the "large-spin" limit. 
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The equations of motion for the Hamiltonians with (J3J) become 
n = p, ± £i jk |B(r,p)| -1 (Bj-a;frfe - Bkufc) , 

Pi = -w?rj ± eyfc |B(r, p)\~ 1 (B j ujfp k - B k uj\p 2 ) . k € {x, y, z} (5) 

This is a non-linear system of six equations, and the search for periodic orbits is 
not easy. We have determined them by a Newton- Raphson iteration, employing 
the stability matrix that enters the amplitudes in the trace formulary Special 
care must be taken at mode-conversion points where B(r, p) = and hence the 
right-hand sides of (jjj) diverge. Some of the periodic orbits, however, can be 
found more easily. This follows from the fact that the three planes r k = Pk = 
(k = x,y,z) in phase space are invariant under the Hamiltonian flow. The 
equations of motion for these two-dimensional orbits are almost harmonic: 

h = Pi T tijk sign(B fe ) u?-rj , 

Pi = -^n T sign( J B fc ) ujfpj . k € {x, y, z} (6) 

Harmonic solutions of @ can be obtained analytically by ignoring the factors 
sign(i?fc) that make the equations nonlinear, and by accepting only solutions 
with the correct constant values of sign(Bk). In total, we find six doubly- 
degenerate planar solutions with different frequencies uf- (i ^ j) given by 




These solutions have the form of ellipses, sketched in the left part of Fig. [| 
which for k — ► shrink to librations along the Ti axes with the original fre- 
quencies Wi. An important fact for these orbits is that B& ^ for all k^O. 

It turns out that the frequencies ojf^ , whose dependence on n is displayed in 
the right-hand part of Fig. for k — 0.1 closely agree with the peak positions 
found in Fig. [|. We cannot resolve all six main peaks; a slightly better resolu- 
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Figure 3: Left: A schematic plot of a pair of the two-dimensional periodic solutions 
of (^|). Right: Frequencies ui^ (Q) of the six planar orbits versus spin-orbit coupling 
parameter k (using the same parameters Wj as in Fig. ti.) 
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Figure 4: Oscillating part Sg (E) of the level density of the 3-dimensional anisotropic 
harmonic oscillator (same parameters as in Fig. but Gaussian averaging range 
O.bhui). Solid lines: quantum-mechanical results; dotted lines: semiclassical results, 
for k = 0.1 using the six planar orbits (first repetitions only) in the trace formula. 



tion is found amongst the peaks corresponding to the second harmonics. Some 
additional small peaks not corresponding to any of the lu^ (e.g., the tiny peak 
at T = 10.1) can be explained by genuine three-dimensional orbits.il 

We now calculate the oscillating part Sg (E) of the level density by in- 
cluding the lowest harmonics of the six planar periodic orbits in the trace 
formula. The ^-correction terms for the planar orbits give a constanL shift 
AS^ = (Xserry + ^nn)T^ = kitcj^/u^ of the actions at all energies£3 The 
result for Sg (E) obtained with k = 0.1 is shown in Fig. || by the dotted lines. 
The quantum-mechanical results are shown by the solid lines. The Gaussian 
averaging range was here 0.5 hw x in order to suppress the contributions from 
second and higher harmonics. For comparison, we also show the curves for 
k = which exhibit a completely different shell sctructure. A good agreement 
is found also in the presence of the spin-orbit interaction, except for low ener- 
gies where semiclassical approximations usually are worst. This shows that the 
planar orbits are the most important ones for the gross-shell structure. The 
numerical calculation of the three-dimensional orbits and their properties is 
not difficult, as long as they are not crossing the mode-conversion surfaces. A 
study of their influence on the finer details of the shell structure is in progress.c3 
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